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ABSTRACT 

Some approaches to 2d gravity developed for the last years are reviewed. They are 
physical (Liouville) gravity, topological theories and matrix models. A special at- 
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1 Introduction 



After Polyakov proposed in 1981 the way to deal with non-critical strings |l[], there were many 
attempts to get some concrete results. However, it still remains one of the most important 
and challenging problems in modern string theory. Indeed, instead of ordinary (string) matter 
theory Polyakov considered the theory of matter fields coupled to two-dimensional gravity. 
This is why the investigation of such theories is of great importance. 

Actually, the revival of the interest to these theories occured after 1988, mainly due to 
the calculation of the critical indices proposed by some groups both in light cone p[ and in 
conformal ^ gauges. These groups have given strong evidence in favor of the existence of 
self-consistent non-critical strings out of "magic" interval of the central charges of the matter 
fields of the theory 1 < c < 25. 

The next serious sucsess was achieved in the fall of 1989 when there were invented some 
new "non- field" viewpoints to the theory of 2d gravity which have given efficient tools to 
deal with the matter theories with the cental charge c < 1. Below I try to review this new 
development as it is observed from nowdays. 

Indeed, in section 2 I just sketchily describe the "field theory" approach of 0, |^, i.e. 
Liouville field theory. But this approach, which I refer to "physical gravity" is still poorly 
developed, so the next sections are devoted to some new frameworks like topological theories 
of different types (section 3) and matrix models (sections 4 and 5). Some concluding remarks 
can be found in section 6. Let us emphasize that we intentedly present (in section 5) only 
sketchy review of matrix models in external fields. This very important approach deserves a 
separate paper. We hope to present it in the second part of the review, which is to be devoted 
to matrix models in applications to wider set of physical problems (like gauge theories). 

In fact, the main part of the review is devoted to matrix model approach, the first, as it 
is the main point of authors's interest and, the second, as just within this approach the most 
striking, powerful and general results were obtained. This review is goaled to non-specialists, 
therefore, almost all tedious technical details are omitted, but the text is provided by a number 
of examples and simple calculations. To get more detailed discussion, one can use references 
to original papers presented in all necessary places. I would like to stress that throughout 
the paper I mainly follow the approach developed by our group. Therefore, the references 
to different approaches are poorly presented. They can be found in two very detailed recent 
reviews 0, |^ . 

2 Physical gravity 

2.1 Liouville theory 

2.1.1 Classical Liouville theory 

In this section we will briefiy consider the Liouville theory having in minds that it describes, 
being properly interpreted, the theory of physical 2d gravity. To begin with, let us consider 
the theory with the canonical Liouville action: 
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Sl = — I dz 



(2.1.1) 



where the field depends on variables z, z, standard complex variables in the Euclidean space- 
time and d = d = This theory has two main features - it possesses conformal 
invariance at the classical level and it has no constant-valued ground state solution, which 
results in disastrous infra-red singularities. 

As for the first property, naively the theory is not conformal, since the trace of the stress 
tensor is equal to Tzz = -^e^^"^. The chiral component of the stress tensor for the action 



pXlD is 



(2.1.2) 



and so for T^z- However, it is not holomorphic, even on the equations of motion, since 



dT,„ 



V2 



(2.1.3) 



This naive stress tensor can be improved by adding the term which does not effect to the 
conservation law: 



Now, taking E = we get 



2^ ^' 72/3 ^ 



(2.1.4) 



(2.1.5) 



and 



dT = ^e^'''^d(P + -^d^d<P = 0, (2.1.6) 

with using the equations of motion. 

In fact, S adds to the action just full derivative, proportional to dd(j), which might be 
essential only when considering non-trivial boundary conditions like the presence of vacuum 
charge |P, 0, H. This term in action can be invariantly written as 

^R<P, (2.1.7) 

where R = dd log p, p = Qzz are the two-dimensional curvature and metric respectively. The 
Poisson brackets of the stress tensor ( p^.l.5| ) produce the Virasoro algebra[| with conformal 
anomaly presented already at the classical level: 

cci = J,. (2.1.8) 



^It is evident from vanishing of trace of the improved stress tensor on the equations of motion. 
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2.1.2 Quantum Liouville theory 



Now let us treat the quantum corrections to the LiouviUe theory. Indeed, the theory is ill- 
defined due to infra-red singularities. Nevertheless, one can try to calculate the quantities 
not sensitive to the infra-red divergences. Say, one can consider ultra-violet renormalizations 
of the action ( |2.1.1|) with the curvature term ( |2.1.7| ). 



Indeed, it is not difficult to take into account all ultra-violet divergent diagramms, as 
in two dimensions these are only tadpoles. The calculation is rather trivial, and results in 
multiplicative renormalization of the coupling constant of the exponential in the action 

([2XTD by the factor A^^/^ . it is the same as shift of the vacuum expectation value of field 
(j) by the quantity 5 = —\^p\ogA. Let us note that this shift depends on ultra-violet cut- 
off parameter A which should be made dimensionless by some infra-red cut-off parameter. 
This parameter can be naturally chosen to be the order of magnitude of the metric p, as 
it is metric that cuts-off the theory at large distances in curved space|. Thus, we can put 
— ^ — -\/2/51ogAp. It just results in the renormalization of the coefficient in the curvature 
term (|2T7D : 



-d4)d(f) -> --(0 + S)dd{(f) + S) = -d4>d4> - ^ddS = -d^dcf) + V2l3R(t). (2.1.9) 
2 2 2 2 



It results in the cental charge equal to 



l + 6(/3+i)2. (2.1.10) 



Now let us justify (or illuminate) this result. Let us suggest that one should still preserve 
the conformal invariance at the quantum level. It means that the exponential term in the 
action should be of the dimension (1,1), i.e. should have canonical dimension. From the other 
hand, the conformal theory with the action ( p.l.l|) and additional aRip term has the central 
charge 0, |^ c = 1 + 3a^, and the "naive" dimension = ^ of the operator : exp{70} : |^, 
which it has in the gaussian model, changes to = \'y{a — 7). As this dimension is to be 
unit, one gets the expression for the coefficient a: 

a = v^(/3+i),/3 = ^±/f~^, (2.1.11) 



which coincides with that in the formula ( |2.1.9| ). It leads to the same expression for the 
central charge ( |2.1.10|) . 



Certainly, this consideration was rather a motivation than the rigid result. Indeed, the last 
reasoning is not too correct, as we treated only chiral objects throughout the discussion, but 
one can separate chiral and anti-chiral parts only on the equations of motion. This approach 
can be done completely correct only in the case of = 2 Liouville theory P|. 



^It again can be turned to the matter of boundary conditions. 

•^Throughout the review we omit the sign of normal ordering where this can not be misleading. 
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2.1.3 Sin-Gordon viewpoint 



Now we will say some words on how to deal with the Liouville theory in a more correct 
way. Let us note that the exact computation of all ultra-violet contribution would lead to the 
complete final answer in the absence of infra-red divergences. Indeed, there are no dimensional 
parameters, excluding UV cut-off parameter, and, therefore, any non UV contributions to the 
effective action are forbidden by the dimension argument. Unfortunately, this is not the case in 
the Liouville theory as it has terrible infra-red divergences, and, therefore, essentially depends 
on the IR cut-off parameter. These IR divergences are due to the absence of finite-density 
constant classical solution of the theory, i.e. due to the absence of classical vaccum in the 
theory. In principle, this ill-defined behaviour should be expected from the very begining as 
the theory is the conformal one, and, therefore, has no good particle content, finite S-matrix 
and etc. It implies that, to deal with the Liouville theory, one should consider this as the 
theory given in a fixed point of some larger theory. The natural choices of this larger theory 
are the sin- or sinh-Gordon theories, depending on whether (3 is real or pure imaginary^. It 
was investigated in the paper |]10|, where the authors worked with the sin-Gordon theory (i.e. 



P = iP is pure imaginary) with the action^: 

SsQ = ^ f d^z \d(t)d(l) + + 7_e-*^'^'^ + laRcj)] . (2.1.12) 

For the sake of convenience, we also replaced a ^ a = —ia in the action, i.e. c = 1 — 3a^. 



The authors of JIOl generalized field theory calculation above to this more complicated case 
and have calculated /^-functions perturbatively in 7-1- and — 1. They have observed that 
q;//5 is a RG invariant and calulated /9-functions. These were: 

^^,& = ia^/3V7-, (2.1.13) 

The theory appears to have both an infrared and an ultraviolet fixed points in the region 
of approximation. There are two essentially different regimes which were realized. One of 
them corresponds to RG flow from the free fleld theory with a background charge a+ to the 
free fleld theory with a different (and larger, in complete agreement with Zamolodchikov's 
c-theorem [TT|) background charge In the other regime, one flows from an UV flxed point 



with 7+ = 7_ = to an IR flxed point, where 7_ = 0, but 7+ 7^ 0, i.e. one of the exponentials 
(marginal operator) survives in the IR limit, leading to the Liouville action. 

From the equation ( |2.1.13| ) one can determine the values of parameters corresponding to 
this flxed point: 

■^The consideration above was equally correct for any complex j3. 

^It is essential for calculations to have two different constants 7+, 7_ in action, though one of them can 
be removed by a shift of the field. This shift is often infinite in the cases under consideration. 
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a 



/a2 



P = — ±J— + 1. 



(2.1.14) 



This is consistent with the connection ( p. 1.11 
The manifest calculation of 



thus justifying our results. 



Toll also reproduced the values of dimensions which were 



obtained above, and the generalized central charge, which on the RG trajectory is given by 
some monotonic decreasing function (in accordance with Zamolodchikov's c-theorem 
the fixed point describing Liouville theory reaches the value ( p.l.lO|) . 



rT|]) at 



2.2 Liouville gravity 

2.2.1 Polyakov-KPZ-DDK approach 

In the previous section we considered Liouville theory with no reference to the main purpose 
of the review, that is, the decription of 2d physical gravity. In the section we are going to 
discuss 2d physical gravity itself. 

One should start from the original covariant Polyakov action for the string in fiat D- 
dimensional euclidean target space: 

S = ^l d'^^g'^''d,X%X\ z = 1, D. (2.2.1) 
The main purpose is to calculate the correlation functions through the following path integral: 

< ViV2...Vn >= J VXVg^^'e-^ (2.2.2) 



with the action ( |2.2.1| ). We will proceed quantizing this system in conformal gauge, though 
there is an equivalent light-cone gauge approach The conformal gauge approach was 
proposed in the original paper by Polyakov and was developed in [^, where the central 
charge and weights of operators were obtained. 

Thus, we fix conformal gauge g^^ = e'^g^y, where g^^ is a background metric. The gauge 
fixing procedure also introduces reparameterization ghosts 6, c, which are almost inessential 
in our consideration. Therefore, the complete action consists of three pieces: the ghost action. 



the Liouville action (|2.1.1| ) with the proper curvature term and the matter field action. As the 



complete action should be conformally invariant and should not depend on the background 
metric g^y, the complete central charge must vanish]^. 

In fact, the ghost stress tensor contributes -26 into the full central charge, and the matter 
fields do D. Certainly, one can consider, instead of (|2.2.1|) , an arbitrary conformal matter 
theory, with an arbitrary central charge instead D. We will do so to keep the most general 
expressions admissing non-integer values of the central charge. 

Now we can calculate the restriction due to zero full central charge, using also the manifest 
expression for the Liouville central charge (|2.1.1(]|) : 



'This means that the theory ought to be topological one, see the next section. 
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c,n + Cgh + CL = ^ c„ + (-26) + (1 + 3^2) = 0, (2.2.3) 

i.e. 



(2.2.4) 

and a was determined in the previous subsection. 

It is evident from the expression ( [2.1.11 ) that the condition < 8 leads to the complex 



values of (3 and the Liouville central charge (|2.1.1CI| ), i.e. to non- unitary theory. It repro- 



duces notorious limitation on the central charge of matter (i.e. on the dimension of the target 
space): 

Cm < 1. (2.2.5) 

In fact, complex (3 just means that the chosen vacuum of the theory is unstable, and this 
vacuum is to pass to a different ground state. It can be understood also from the tachyon 
mass square which is equal to This means that the tachyon is absent and the vacuum 

state is stable only provided by the same condition < 1- 

There is another limit when cl is also real. It corresponds to < 0, i.e. c > 25. This 
implies that (3 is pure imaginary, or, what is the same after redefiniton (p —>■ the wrong 
sign of the kinetic Liouville term. It seems to spoil the unitarity of the theory too. 

But there is a special case of Cm = 25, when one just reproduces the critical string with 
25 space coordinates (matter fields) and one time coordinate with reversed sign in the kinetic 
term (Liouville field). Therefore, in such a way, we recover the mechanism to get Minkowski 
space from Euclidean theory, which is applicable for the critical string! Simultaneously it 
means that the "physical" dimension of non-critical string is Cm + 1- 



2.2.2 Dressing of operators 

Now we are going to discuss operator content of 2d gravity. As the first step, let us note that 
the functional integral in gravity theory should be independent both of background metric 
gfj,i, and of physical metric g^u = g^uC^^^^ 

The first property implies zero full central charge of the theory and was already used to 
connect Liouville exponential with the matter central charge Cm- The second requirement was 



used in the previous subsection to calculate central charge of the Liouville theory (|2.1.10|) . 

Now we will consider matter field in the same context. These fields should be also renor- 
malized (dressed) by the Lioville field due to the interaction with gravity. To calculate this 
dressing, we require for correlation functions of the operators to be independent of metric, 
instead of doing very tedious perturbative calculations. 



^This renormalization of the exponential factor - or the Liouville field - is a direct consequence of our 
requirement for this exponential to have the dimension (1,1). The same phenomenon can be described as 
well by saying that the scalar Liouville field changes the trasformation law and is not scalar after taking into 
account the conformal anomaly (see more detailed discussion in M). 
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Thus, let us consider a matter operator V which has a dimension A. It is used to suggest 
that, as a result of gravity interaction, it is dressed by exponential of the Liouville field: 



Vdr = e'^V. (2.2.6) 

Actually, we would have to consider more general operators which include Liouville secon- 
daries (and maybe even ghost contributions). Fortunately, it is not necessary to do in the 2d 
gravity case. This is due to huge invariance of the theory. The drastic reduction of the space 
of vertex operators is a general phenomenon occurring whenever a string model is build from 
CFT. Let us remind first the situation with the critical string. In the simplest treatment, one 
neglects the Liouville field at all, but this is not the only thing one does. There are three 
additional requirements to physical vertex operators in critical string models: 

a) They do not contain ghost fields (in certain and natural picture)^. 

b) They are Virasoro primaries (with respect to "full" Virasoro). (All the descendants, 
present in CFT are "gauged out" . The reason is that coupling to 2c? gravity implies gauging the 
Virasoro algebra and thus all the descendants are eliminated as gauge-non-invariant operators. 
They really decouple from the correlators^) . 

c) They are integrals (in the picture consistent with a)) of operators of conformal dimension 
one0- 

Naive generalization of these three principles to the case of non-critical string has been 
suggested in p, with the only modification made at the point b), where one is supposed 
to take primaries of the full {matter + Liouville) Virasoro algebra, in the form of {matter 
primaries) x {Liouville primaries) and Liouville primaries being pure exponentials. It is, how- 
ever, not quite true that such a simple generalization is valid. A problem arises at the very 
beginning - at point a). In order to derive honestly this requirement one should prove that 
in every BRST cohomology class (i.e. for every physical state) there is a ghost free repre- 



sentative. According to [jT2|, |13| this is not true in non-critical case: there are two additional 
representatives of the cohomology classes, which unavoidably (and non-trivially) contain ghost 
fields. In the particular case of c = 1 model, the first ones were interpreted in [2] as a ground 
ring, another one (the "ghost number two" ) still has no nice interpretation. It seems, however, 
that the subsector, described by principles a),b),c) is closed by itself under OPE (modulo de- 
scendants, i.e. fields vanishing in the correlation functions). Another delicate point is that the 
requirement c), i.e. that the full dimension is one, permits two different choices of Liouville 
primaries associated with a given matter operator. In what follows we consider a subsector 
with one specific choice of these two. This subsector also seems to be closed under OPE in 
the above sense. 

Thus, we restrict ourselves to this subsector and consider only the operators of the form 
( 2.2.6|) , more general discussion can be found in [14|. Let us calculate the value of the exponent 



7 in ( ^.2.61 ) requiring (1,1) dimension of the dressed operator (it was already explained in 
different words that it just expresses the possibility to integrate it over 2d surface with no 
breaking conformal invariance). Then we obtain: 



^Except for one delicate point, concerning dilaton operator. 
^It is just no-ghost theorem. 

^°Let us note that other (equivalent) pictures may be obtained by multiplication of dimension one integrands 
(constrained by a) and b)) with ghost field instead of integration. 
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A + -7[a - 7] = 1, 



(2.2.7) 



I.e. 



1 



'y = -[a±Va^ - 8 + 8A). (2.2.8) 
Now let us say some words on the choice of sign in this expression (and that in the Liou- 



ville exponent, see ( 2.1.11 )). Hereafter, our prescription will be to choose because of 
arguments proposed in There was argued that this choice selects out the local operators 
corresponding to the states with positive energy. We will comment more this point a bit later. 

In the next subsection we use this expression to determine scaling properties of dressed 
operators. 



2.2.3 String susceptibility 

As correlation functions of physical operators in 2d gravity do not depend on metric, as well 
as on the positions of operators (since the integration over locations of all operators is to be 
presented), all correlation functions are just numbers. We will return to the question in the 
next section and now let us only note that we still need some characteristic like dimension, 
which differs among operators (and correlation functions). The simplest way to introduce such 
a characteristic is to fix the area of the surface A and investigate A-dependence of correlation 
functions. It is clear that this characteristic is of great importance, as it is that which governs 
the grow (or shrinking) of the surfaces and controls their stability. It is the general feature of 
all string theories that their partition functions exponentially decrease with area. The index 
which really differs among string theories and governs the behaviour of their ground states is 
pre-exponential factor. 

We will now calulate this factor. Let us consider the path integral for 2d gravity (=non- 
critical string) theory with fixed area A: 

Z{A) = J P(/.P(matter)P(ghosts)e-^M (^J d^^y^e^^"^ - A^ , (2.2.9) 

where we included matter and ghost action terms into the measures and took into account 
that the area of the surface is nothing but the integral of the physical metric g^j^jy — g^v^ 
over the surface. 



Though the integral (|2.2.9|) is very complicated, its A-dependence can be easily determined. 



Indeed, let us do the substitution — > + log A. Then, the exponential in the Liouville 

action turns into the exponential term g-constxA ^j^g functional integral (due to the 5- 
function in ( p.2.9|) ). Therefore, the only contribution of the action to the pre-factor is from 
the curvature term: 

^ J d^^fgaR<P ^^J d'^fgaR<P + log A (i- / d'^fgR^ . (2.2.10) 
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Measure in the path integral ( |2.2.9| ) is evidently invariant with respect to constant shifts. 
Then, using the identity 6{Ax) = A~^6{x) and the Gauss-Bonnet theorem ^ / d?^\/gR = 
I — g, where g is the genus of the surface, one finally obtain^: 

Z{A) ~ ^-;%(^-^')-^e-constxA _ ^{7,t,-2)(i-g)-ig-constxA^ (2.2.11) 

i.e. 

Istr = 2--^ = ^(^Cm-l- v^(c^-25)(c^-l)) . (2.2.12) 
This index jstr is called string susceptibihtyQ. 



2.2.4 Weights of operators 



In the previous subsection we have calculated the index describing the dependence of the 
partition function on area A. Analogous dependence of correlation functions can be also 
investigated, and defines the index which is a proper characteristic of operator in 2d gravity. 
More concretely, the index, which is called the weight of operator, can be determined from 
the normalized one-point correlation function: 



< V >A= / r'0r'(matter)P(ghosts)x 

(2.2.13) 



xe 



l-h 



Physically, this index controls the behaviour (growth) of string surface when the operator is 
inserted. Therefore, the standard notions of relevant, marginal or irrelevant operators still 
make sense. In particular, relevant operators correspond to /i < 1 and dominate in the IR 
(large A) limit. 

To calculate h, we do the same steps as in course of calculation of the string susceptibility 
and use the manifest expression for the dressing exponent ( [2.2.8| ): 

< V >A , = AV-^^, (2.2.14) 

A 

i.e. 



v2/5 v25 — Cm — V 1 — Cm 

Thus, the classification of operators into relevant, marginal and irrelevant ones does not 
change for dressed operators. Indeed, the conditions A = 1, A > 1 and A < 1 immediately 
imply the same conditions for the dressed weights h in ( ^.2.15] ). 

^^We do not concretize the constant in the exponential because of its non-universal nature, in contrast to 
the index jstr- 

^^Let us remark that jstr has a complex value when 1 < Cm < 25, pointing once again to "the magic 
interval" of values of Cm- 
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Let us note that the expression ( |2.2.15| ) is real and positive for < 1 provided the choice 
of sign is "-" in the expressions (|2.2.8| ) and (|2.1.11|) . This is one more argument in favor of 
that choice. 

Now, to have an exphcit example, let us write down the values of jstr and h for the 
operators in minimal models [l^, |^, §]. The central charge of (p, g)-minimal model is equal 

c=l-^^^, (2.2.16) 
pq 



I.e. 



2b -g| (ooT7\ 

lstr = \ , f- {Z.Z.K) 

p + q-\p - q\ 
The (Kac) spectrum of operator dimensions is given by 

A„,„= 'P'"-^"f-<P-'''' . I<.n<,-1, l<n<p-l, (2.2.18) 
Apq 

i.e. 

_ li»"-H-|p-<;l 2 

p + q-\p- q\ 

This expression should be compared with the indices calculated in matrix models (the below). 



3 Topological theories 



We could conclude from the discussion of the previous section that correlation functions in 
the theory of 2d gravity have a very simple structure. Indeed, if we do not fix the area of 
the surface, then all correlation functions are just numbers, provided by vanishing of the full 
central charge of the theory. This means that it is a topological theory. 

Indeed, there no variables which correlation functions could depend on, since integration 
over metrics is done. Certainly, it implies that the theory is in unbroken gravity phase. 

Let us look at the same argument from a different viewpoint. The routine way to do 
(generally, non-critical) string integral was usually to fix the conformal factor and then to 
remove two remaining components of metric tensor by making use of the general covariance of 
the string action (i.e. to integrate out the diffeomorphism group). After doing this one remains 
with the integral over moduli of the conformal (=complex) structures and with the integral 
over conformal (Liouville) factor (in the case of non-critical strings). It is approximately the 
procedure which was implied in the previous section. 

Now let us change the order of the integrations. Namely, one can integrate over metrics 
from the very beginning. Then, we have an integral (correlation function) which can depend 
only on the locations of the points (as any metric dependence is already integrated out) and 
still possesses the general covariance (which is again correct only in unbroken phase). It can 
be evidently nothing but a number. 
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Thus, we are led to consider the theory of 2d gravity as a topological theory. Certainly, 
there is plenty of different topological theories, but it is naturally to consider only 2d theories. 
Indeed, there are only two sucsessful attempts to formulate a topological theory appropriate 
for the description of 2d gravity. 



The first attempt was due to Witten jl^, who proposed the theory which preserves the 
memory of the moduli space treating as correlation functions some integrals over moduli 
space. More concretely, let us consider the holomorphic linear bundle J-'i with the fiber being 
the cotangent space (complex one dimensional vector space) to the Riemann surface S at the 
marked point Xj. Then we can take as correlation functions the intersection numbers of these 
bundles. It means that we should integrate the first Chern classes of the exterior products 
of these bundles over (Deligne-Mumford) compactification A4g,n of the moduli space of the 
curve S of the genus g (see also fl^ 



<a,,...a,„>^ f nci(-^.)'^'^ (3.1) 

JMg,n j 

These correlation functions are indeed numbers (positive rational) which vanish unless Qg — 
Q + 2n = 2Yl,di. They are related by a number of recursion relations which are equivalent to 
the Virasoro constraints of the Hermitean one- matrix model in the continuum limit |19|. This 
allows one to identify the two theories. Unfortunately, it is not so immediate to generalize 
the construction to multi-matrix model case, and the work in this direction is not finished yet 
(see, however, pT|). 



Another sucsessful construction of the appropriate topological theory was based on the 
follwoing general idea. Let us consider the theory whose stress tensor (or action) is anticom- 
mutator of something with a BRST operator Q: T = {*,(5}, (5^ = 0. Then, this stress tensor 
describes a topological theory on the subspace of the physical states (i.e. those defined by the 
condition Q\phys >= 0): 

T\phys >= 0. (3.2) 

In fact, since coordinate variations are generated by the stress tensor, any variation of a 
correlator is propostional to the the same correlator with insertion of T. Therefore, correlator 
of the physical fields, due to the condition (|3.2| ), does not depend on the locations of the 
points and is a number. 

A concrete realization of this general idea has been proposed in |^ (see also [^), where 
twisted N = 2 superconformal theory was investigated. It was formulated in abstract terms of 
OPE, with no any Lagrangians. The Lagrangian formulation, namely, N = 2 supersymmetric 
Landau-Ginzburg theory has been proposed in . 



The advantage of the theories of described type is that they can be investigated in some 
details, and they are applicable to the description of plenty of matrix models, in contrast to 
the theories of Witten's type. Because of the lack of space we can not develop here all these 



ideas, and refer to the reviews and papers containing all necessary references p5|, p6|, p7|. 



4 Discrete matrix models - the first step 
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4.1 A role of matrix models 



Now let us say that the resultes mentioned in the section 2 are the only ones obtained in the 
framework of physical gravity. But in the Fall of 1989 a drastical progress in 2d gravity was 
achieved in some other, rather uncommon for the field theorists, frameworks. In particular, an 
old idea proposed in to replace the path integral of 2d gravity over metrics ( p.2.2|) by the 
sum over all possible triangulations of surfaces (i.e. reformulate it in target space terms) have 
led, after all, to what is called nowdays matrix models. Really it consists of three different 
components - double scaling limit p9|, constraint algebras satisfied by the partition 

functions of matrix models |1^, ^ and integrability properties |^ . 



Up to now, the approach relies only on the original "physical" arguments of triangulationQ 
with sucsessive taking the (double scaling) continuum limit and on the comparison of the very 
poor predictions obtained from the physical gravity with those within matrix model treatment. 
These predictions include weights of operators in the theory and string susceptebility, as well 
as the calculations of 3-point functions (and some n-point functions in very simple channels). 
In fact, it means that the comparison is at the level of the spectra of the theories. It implies, 
more or less, that the matter of quantum measure is still out of the scope of the present 
investigations in physical gravity. In fact, though there are some distinguished measures in 
the gravity path integral (induced by the free field representation, or Polyakov's one, or ...), 
it is still necessary to fix somehow these as well as the normalization factor depending on 
the genus of the (Riemann surface) string world sheet (i.e. on the order of the perturbation 
theory) . 

At this point, one can reverse the line of reasoning and consider matrix models as a 
definition of these entries. To be a good definition, such a theory should be defined naturally 
over all genera simultaneously. Indeed, I will try to demonstrate below that matrix models 
have to do with the infinite Grassmannian, which, in its turn, naturally describes the Universal 
Moduli Space (UMS, the space of all Riemann surfaces of all genera). Therefore, matrix 
models do really give the example of the theory with required properties. 

Indeed, we will see that the partition function of a matrix model is a r-function and 
corresponds to a concrete point of the Grassmannian (corresponding to a Riemann surface 
of infinite genus) which implies that after doing the path integral (|2.2.2| ) over UMS, one 
reproduces the answer corresponding to a point of the same UMS (see also [Bl, BR PB|). It 



resembles the phenomenon realized some years ago in 2d quantum integrable systems, 
when the exact quantum correlators satisfy classical integrable equations^. 

Thus, it is expected that matrix models give a correct description of the non-critical 
strings with c < 1, therefore, better understanding of their properties is necessary. I will try 
to describe briefiy these properties in the next two sections. Indeed, now we can do even 
more - we establish the connection of matrix models with some other recent approaches to 
2d gravity which can be also solved up to the very end (in contrast to physical gravity). It is 
possible to reach the proof of complete identity between all these theories. This should not be 
very surprising as all these theories calculate, in essence, the same cohomologies and refiect 

^•^which are so approximate that even can not give any physical interpretation of multicritical points 
^^The analogy is even closer as in string theory one should integrate some objects having r-function sense 
giving rise to matrix model partition function which is a r-function again, and in the authors obtained 
the classical integrable equations for a functional of correlation functions which was also a r-function. 
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different general properties of the same final answer. 

Now we just start from the matrix integral, the standard arguments in favor of it can be 
found in plenty of reviews @, H, ^ . Let us note that all the specific properties of matrix 
models will be demonstrated in this section for the example of Hermitian one matrix model. 
More complicated examples (in this and the next sections) will be labelled by asterisque and 
can be omitted without disturbing the main line of reasoning. I omit many subtle points and 
all evident (but maybe technically interesting) generalizations of the considered examples, 
refering to the original papers. Many of these generalizations can be also found in two recent 
reviews 13^, EO 



4.2 Integrability in matrix models 
4.2.1 Toda chain equation 

Let us consider the integral over Hermitian n x n matrix with arbitrary potential: 

Zn{tk) = {Voluin)n\}-' J[dH]exp{-TTV{H)}, V{H) = -Y,tkH\ (4.2.1) 

4 = 

where [dH] is the proper Haar measure. The standard procedure of doing this integral is, 
indeed, to fix some simple potential (polynomial of a finite degree), to write down the equations 



satisfied by the partition function (usually there are few equations as the integral ( [4.2.1 



depends on some coefficients in the polynomial potential) and to take the double scaling limit 



in these equations. Instead, we consider the partition function ( |4.2.1| ) as the function of all 
coefficients t^s and demonstarte that all the crucial properties of matrix models (integrability 
and constraint algebra) can be exposed before taking the continuum limit. Moreover, in a 
sense, these properties are easier observed in the case. 

Thus, to begin with, let us demonstrate that the partition function ( [4.2. 1|) is nothing but 
the r-function of the Toda chain integrable hierarchy [^, We can rewrite ( [1.2. 1[ ) as the 
integral only over eigenvalues with angular variables integrated over: 

Mtk) = {n\)-' Jl[dh,A\h)exp{-J2t^h1}. (4.2.2) 

i i,k 

Now one can apply the standard machinery of the orthogonal polynomials To do this, 
let us define the polynomials by the orthogonality condition 

m^^p,.p,) = |p.wPAl^^'^''ih=s,e^^'^ (4.2.3) 

where e'^*'^*^ are the norms to be determined, the polynomial normalization being fixed by the 
unit coefficient in the leading term: 

Pi{h) = J2aijh\ aii = l. (4.2.4) 
Using these polynomials, one can rewrite ( [1.2.2[ ) as 
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n-1 



Zn = (nl)-' Jlldh,detPk-i{h^)detPi.,{hJexp{-J2t,h1} = H e'^^^*). (4.2.5) 

i i,k 1=0 

Now we are going to prove that this partition function is nothing but a r-function of Toda 



chain hierarchy. In fact, we derive only the first equation, refering to |^ for the general proof 
and any details. 

To begin with, let us consider the scalar product < nlhlm > for different m. It is evidently 
zero whenever m — n 0, ±1 (e.g., < n\h\n — 2 >=< n\n — 1 > + Y^k>i < n\n — k >= 
by the conditions and ( ^^ ). It means that 



hPn{h) = Pn+l{h) - PnPn{h) + RnPn-l{h), (4.2.6) 

where the coefficients p„ and Rn should be determined. The latter can be trivially calculated 
using ( |4.2.tj ) twice: < n — l\h\n >=< hPn-i,Pn >=< n\n > +J2k>o < "-1^ ~ ^ >=< 
n\n >=< Pn-i, hPn >= Rn < n — l\n — 1 >, i.e. 



i?„(t) = e'^"-'^"-!. (4.2.7) 
Indeed, the recurrent relation ( |4.2.6| ) already hints at the Toda chain, as, introducing difference 



L-operator with the matrix elements 



^m,n—l Pn^m,n ~l~ Rn+l^m,n+l 



(4.2.8) 



with eigenfunction (Baker-Akhiezer function) "^nih) = exp{ — ^V{h)}P]\f{h), one can trivially 
see in ( |4.2.6| ) the action of Lax operator of the Toda chain [^, where h is, as usual, the 
spectral parameter: 



L*„(^) = h^n{h). (4.2.9) 

Now to get the equations of Toda chain it is sufficient to obtain the correct action of the 
operator of the derivative with respect to the first time [^. But more illuminating way is 
to derive the Toda chain equations [^]. As an example, we demonstrate here the derivation 
only of the first equation in the hierarchy. To do this, let us differentiate < n\n > with respect 
to the first time using (^X^ ), (|g and (^^■. 



d /dP \ 

— < n\n >= 2 ( — ^, P„ ) - < n\h\n >- 



I.e. Pn 



(4.2.10) 

where dot means the derivative with respect to the first time. Now we differentiate the same 
quantity once more and, using (^4.2.3|) , (|4.2.4| ), (|4.2.6| ), (|4.2.71 ) and (|4.2.10| ) and integrating by 
parts, we obtain finally the well-known Toda chain equation: 



(9^ / dP ^ 

^ne"" + {^nfe^" = — <n\n >=< n\h^\n > -2 (-^, hP^ 

= (pI + Rn+i + Rn) e^" - 2Rne^\ 



(4.2.11) 
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i.e. 



(pn = Rn+l -Rn = 6'^"+^"'^" - e'^""'^"-!. (4.2.12) 

This equation can be rewritten in the Hirota bihnear form leading to the notion of r-function, 
key object in the theory of integrable hierarchies. In fact, let us introduce this by the formula 

e^"W = (4.2.13) 



Then, one can rewrite ( [4.2. 12[ ) in the Hirota bilinear form as 



rn{t)^Jn{t) - = r„^.i(t)r„_i(t). (4.2.14) 

Indeed, the main sense of the notion of r-function is the statement of existance the unique 
function depending on infinite set of times and satisfying the infinite hierarchy of equations. 
What we demonstrated was the only first equation of whole hierarchy, see more exhausting 



consideration in E3, WL 



Hereafter, we will not address to any equations and discuss only r-functions. What we 
are going to show for the next two sections is that any matrix model partition function is a 
T -function. Indeed, using ( |4.2.13| ), ( |4.2.5| ) and putting 



To = 1, (4.2.15) 

one can obtain for one matrix model: 

Zn = - = r„. (4.2.16) 

To 

4.2.2 Determinant representation 

Certainly, to demonstrate that the partition function is a r-function one should not check 
all infinite hierarchy of equations. During our consideration in the next sections, we obtain 
some effective and convenient representations for r-functions. Now let me just prove that the 
r-function of the Hermitian one matrix model can be presented in the specific determinant 
form whichj defines, actually, Toda chain r-function. 

Let us rewrite the orthogonality condition ([4.2.3|) in the "matrix form". Namely, we 
introduce matrix A with matrix elements a^n determined in ( |4.2.4| ), so-called moment matrix 
C whose matrix elements are determined by 

C^, = ( dhh'+^-^e-^^^^ (4.2.17) 



and the diagonal matrix J whose diagonal entries are just e'''". Then one can rewrite ( [4.2.3|) 
as matrix relation^ ACA^ = J {A^ means transponed matrix), and taking the determinant 



of both sides of the relation, one obtain 
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This equation is nothing but Riemann-Hilbert (factorization) problem - see |44, ^ for details. 
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Zn — det Cij 



(4.2.18) 



the conditions 



and 



dtu 



dtf 



c — c 



(4.2.19) 
(4.2.20) 



C„ = d'^-^Cii = (4.2.21) 

being satisfied (as it is evident from tlie manifest form of the moment matrix ([4.2.17|) ). Thus, 
one can write finally 



Z„ = det d'+^-^C. 



(4.2.22) 



Note that the conditions ( [4.2.19| ) and ( |4.2.21| ) is to be satisfied for any (Toda) KP hierarchy, 
while ( [4.2.201) determines Toda chain hierarchy [^. We will return to the question later. 



4.2.3 * Multi-matrix models 



In this section we discuss the integrability in multi-matrix models considering at the beginning 
the standard multi- matrix model considered in |4^. Namely, its partition function is given 
by the matrix integral]^ 



/ nM^.]exp|-Tr 



x:vi(Ho+x:HiH,- 

i=l i=l 



(4.2.23) 



This integral again can be integrated over angular variables, this time using Itzykson-Zuber 



formula 47 



and one obtains 



^trXULW 



(4.2.24) 



p n 



YlYldh 

.1=1 i=l 



{1} 



n dh^^^dk 



.1=1 



p-i 



A(/^ii})A(/^w) exp <! - 5: mhf^) 

i,l 1=1 i 



(4.2.25) 



'For the sake of brevity, we omit all proportionality factors, which can be found, say, in |3£ 
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To apply the orthogonal polynomials in this case, one needs two sets of polynomials Qn{x) 
and Tm{y)i which satisfy the orthogonality relation: 

j dxdyPnix)Qmiy)A{x,y) = 6m,nhmitk,h\ck), (4.2.26) 

where times tk and 4 parameterize potentials Vi and Vp accordingly, while coefficients Ck 
parameterize other potentials. Then, in complete analogy with one matrix model case, one 
can obtain the determinant representation: 

= det C^^ (4.2.27) 



nxn 



where 

Clf = J dxdyx'-^y'-^A{x,y), Cjf^ = C^p> (4.2.28) 



and 



Zn = det d'^^d^-^C^p} (4.2.29) 



as the entries in the determinant satisfy the conditions analogous to ([4.2. 19|) and ( [4.2.21 



which imply that the partition function ([4.2.27[) is the r-function of general Toda lattice hi- 
erarchy Hl[, This property is understood already from the fact that the partition 



function ( 4.2.29 ) depends essentially on two different sets of times (which correspond to posi- 



tive and negative times in Toda lattice hierarchy), with the other times (potentials) in ( [1.2.23[ ) 



just parameterizing the r-function, i.e. the point of the Grassmannian. Certainly, it would be 
more natural to have all times at the equal footing. In fact, we considered up to now only one 
possible generalization to multi-matrix models. But there is one more, from some viewpoints 
more attractive, extension. This is so-called "Conformal Multi-Matrix Models" (CMMM) 



48[ , [49[| , which will be discussed in the next subsection. These CMMM can be described by 
multi-component hierarchies of SL{p + 1) AKNS type and depend on many sets of times. 

For the sake of simplicity, we write down here only two matrix model case. Then, the 
model is given by the partition function (in eigenvalue variables) 



2'{CMMM,2} ^ f 
"i,"2 ~ J 



2 n. 



n n dhii' 
i=iii=i 



n A^(/.«) uihi'^ - hf) exp <! - 5: v.ik 



{On 



(4.2.30) 



1=1 



We will discuss later why this quite specific partition function is especially interesting and 
discuss its integrable structure in more details. Now it is just an example of another multi- 
matrix model which has a rather simple determinant representation 
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det 



d 



C 

dC 



2-{CMMM,2} 

ni,n2 



c 



d' 



,2ni+n2-2 



C d 



C 

Be 

irti +n2-l 



c 



. d 



Qn2-\Q 
2ni+n2—2 



c 



where 



C{t) = / dz exp[-Vi{z)], C(t) = / dz exp[-V2iz)]. 



(4.2.31) 



(4.2.32) 

This representation gives rise to the r-function of (generahzed) SL{3) AKNS system [0, |48[] . 
From this view point, one matrix model corresponds to usual SL{2) AKNS hierarchy, which 
is equivalent to Toda chain |5^. Certainly, analogous equivalence is absent for higher SL{N) 
hierarchies. 



4.3 Virasoro constraints 

4.3.1 Ward identities in matrix models 

In the previous subsection we briefly described the integrable properties of the discrete matrix 
models. In this subsection we pay attention to equally important and general property, 
constraint algebra which is imposed on matrix models. The role of these constraint algebra 
is to fix the concrete r-function. In fact, the integrable hierarchy does not fix the r-function 
uniquely. Moreover, there is a huge space of the solution to the equations of hierarchy, which 
are parameterized by (subspaces) of infinite dimensional Grassmannian. It turns out that 
constraint algebra imposed on r-function fixes it usually uniquely. It means that in order to 
describe matrix model parittion function in invariant way, one should determine what kind of 
r-function is the corresponding partition function (what hierarchy it is described) and what 
constraint algebra it satisfies. Again, let us start with the simplest example of Hermitian one 
matrix model. 



We demonstrate now that the partition function ( [4.2. 1|) satisfies Virasoro algebra con- 



straints (more precisely, its Borel subalgebra). Indeed, these are nothing but Ward iden- 



tities of the model ( |4.2.1|) [p] , p^ , p3| , p4[ . To see this, let us shift the variable H by 



small quantity eW^^^. Then, as it is nothing but the change of variables, the integral 



( |4.2.1|) will be unchanged. From the other hand, this shift results into changing of poten- 



tial E tkH'' EkiH + eHi+^f = ELi tkH'' + T.k>qitk + {k- q)tk-,)H''). This means that 
these shifts are given by the action of operators L'^J; = Yl, ktk p,P producing Virasoro algebra. 

Now we should also take into account the changing of the measure under the transfor- 
mation above. The Jacobian of variation of Haar measure on Hermitian matrices is equal to 
det ll^^^^^Jg^^ll ~ 1 + eTr^U^ = 1 + ELo Tri^'^Tr if As any degree k of matrix H 
can be obtained by differentiating of the exponential of potential with respect to t^, one can 
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reproduce the changing of measure by the operators Y.k=o 'at^t — k' ^^^^ according to 



the definition (|]2]T]) 



—Z -nZ 

dto 



Therefore, the final result of the shift is: 

Zn{t) = Zn{t) Iwith shifted H = (1 + Lg)Zn{t), (4.3.1^ 



I.e. 



L,Zn = 0, q>-l, (4.3.2) 

where operators 

satisfy the commutation relations of Virasoro algebra 

[Ln, Lm] = (n - m)Ln+m- (4.3.4) 

As we consider only Borel subalgebra of whole Virasoro algebra, the question of the central 
charge is out of our consideration. 

Thus, we proved that the partition function of Hermitian one matrix model satisfies Vira- 
soro algebra constraints. Unfortunately, it is not so simple to find out an analogous symmetry 
in multi- matrix models ( [4.2.231 ). So, we are going to suggest some more general way to con- 
struct matrix models possessing given symmetry, which allows one to write down multi-matrix 
models with W-symmetry, and, as a by-product, explains why conformal algebra arises in ma- 
trix models. 



4.3.2 Conformal matrix models 



In fact, now we just demonstrate that the partition function of ( [4.2. 1| ) can be rewritten in 



terms of correlators in gaussian c = 1 confomal field theory and the Virasoro generators ( [4.3.3[ ) 



actually have the well-known form of the Virasoro operators in the theory of one free scalar 

field. Indeed, let us consider holomorphic components of the scalar field 

0(z) =Q + P\ogz+J2^z-' 

[Jn, Jm] = n6n+m,0, [Q , P] = ^ (4.3.5) 

and define the vacuum states 

Jk\0) =0, (r2|J_fc = 0, k>0 

p\0) =0, {N\P = n{n\. (4.3.6) 
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' Half of the stress-tensor[^ components 



^ ^ k>0 



obviously vanishes the S'L(2)-invariant vacuum 

T,|0) = 0, q>-l. 

Then, we define the Hamiltonian by 



a + b=q 
a,b>0 



(4.3.7) 



(4.3.8) 



fc>0 



V{z) 



A;>0 



V2 



d4>{z) 



(4.3.9) 



Now let us solve the Virasoro constraints ( 4.3.2 ) in the general form. Put differently, we derive 
the general parition function possessing the Virasoro symmetry ([4.3.2|) . One can easily con- 
struct a "conformal field theory" solution of ( [4.3.2|) in two steps. The basic "transformation" 



L„(n|e^W . . . 



in 



e^WT„ . . . 



(4.3.10) 



can be checked explicitly and show how the Virasoro generators ( |4.3.3| ) transforms to those 
in gaussian model. As an immediate consequence, any correlator of the form 



(nle^WClO) 



(4.3.11) 



(n counts the number of zero modes, "included" in G - that is the role of the size of matrix 



in ( [4.2. 1| )) gives a solution to ( [4. 3. 2] ) provided by 

[T„G] = 0, q>-l. 



(4.3.12) 



The conformal solution to ( [4.3. 12| ) (and therefore to ( [1.3. 1[ )) immediately comes from the 
basic properties of 2d conformal algebra. Indeed, any solution to 



[T{z),G] = 



(4.3.13) 



is a solution to ( [4.3. 12| ), and it is well-known that the solution to ( [4.3. 13| ) is (by definition of 
the chiral algebra) a function of screening charges in the free scalar field theory given by 



Q± = fJ±-- 

With a selection rule on zero mode it gives[^ 

G = exp Qh 



nl 



(4.3.14) 



(4.3.15) 



^^For the sake of brevity, we omit the sign of normal ordering in the evident places, say, in the expression 
for T and W in terms of free fields. 

^^Of course, the general case might be G ^ Q^^™Q"^ but the special prescription for integration contours, 
proposed in | |55| , implies that the dependence of m can be irrelevant and one can just put m = 0. In fact, the 
problem with the choice of integration contour is a rather subtle point and discussed in |55| . 
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In this case the solution 

= (n|e^(*)expQ+|0) (4.3.16) 

after computation of the free theory correlator, analytic continuation of the integration con- 
tour gives the result 

/n 

^=1 (4.3.17) 

[n 



!Vol U{n))-^ J [dH] exp (- Trt^i/' 



in the form of the matrix integral ( [4.2.1 ). This is the simplest example of CMMM mentioned 



above which already explains the name. 

Thus, we derive Virasoro invariant matrix model of the general type. Certainly, this 
procedure is immediately generalized to other symmetries. 

4.3.3 * Ward identities in multi-matrix models 

Now we are going to derive multi-matrix models invariant with respect to Ty-algebra and 
discuss what is the symmetry of the standard multi-matrix models ([4.2.23|) . 



Let us consider the immediate generalization of proposed procedure. Now we can use 
powerful tools of 2d conformal field theories, where it is well known how to generalize almost all 
the steps of above construction: first, instead of looking for a solution to Virasoro constraints 
one can impose extended Virasoro or W^-constraints on the partition function. In such case 
one would get Hamiltonians in terms of multi-scalai field theory, and the second step is 
generalized directly using screening charges for IV-algebras. The general scheme of solving 
discrete l^-constraints looks as follows: 

(i) Consider Hamiltonian as a linear combination of the Cartan currents of a level one 
Kac- Moody algebra Q 

i7(tW,...,t(-'^'^^))= Y: i'^f^xJk, (4.3.18) 

A,fc>0 

where {Mj} are basis vectors in Cartan hyperplane, which for SL{p + 1) case are chosen to 
satisfy 

I p+i 

MrM, = 5i,-— -, Ea^. = 0- 
^ j=i 

(ii) The action of differential operators W^""^ with respect to times {t^k^} can be now 
defined from the relation 

iy/"\n|e^«*» . . . = (n|e^"*»wS"^ . . , a = 2,...,p+l; t > 1 - a, (4.3.19) 

where 

(4.3.20) 

w(«)(z) = E[MA50wr + --- 

A 
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are spin-a W-generators of Wp+i-algebra written in terms of rank ^-component scalar fields 



56|. 

(iii) The conformal solution to the discrete VT-constraints arises in the form 

^{CAfMM,p+i}^^^ = (n|e^«*»G{g(")}|0) (4.3.21) 

where G is again an exponential function of screenings of level one Kac-Moody algebra (see 
57| and references therein) 

Q{a) ^ Ijia) ^ fcxcjy ^4 3 22) 



{ck} being roots of finite-dimensional simply laced Lie algebra Q. The correlator ( [4.3.21|) is 
still a free-field correlator and the computation gives it again in a multiple integral form 



nci^-^exp - 4'^(/x.a)(>)^'= 

i=l V A,fc>0 



X 

(4.3.23) 



(a,/3) i=li=l 

The only difference with the one-matrix case ( |4.3.17| )) is that the expressions ( |4.3.23| ) have 
rather complicated representation in terms of multi-matrix integrals, the following objects 
will necessarily appear 



n n(4"^ - [det{i7(") ® / - / ® , (4.3.24) 



3 

However, this is still a model with a chain of matrices and with closest neighbour interactions 
only (in the case of SL{j) +1)). 

Actually, it can be shown that CMMM, defined by ( |4.3.21| ) as a solution to the W- 
constraints has a very rich integrable structure and possesses a natural continuum limit 
49| . To pay for these advantages one should accept a slightly less elegant matrix integral with 
the entries like ( [4.3.241) . 

The first non-trivial example is the p = 3 solution to W^s-algebra: an alternative to the 
conventional 2-matrix model. In this case one has six screening charges Q(^"») {i = 1, 2, 3) 
which commute with 

W(2) (z) = T{z) = - [d(f){z)f (4.3.25) 

and 

W^'\z) = j2{fi,dcP{z))', (4.3.26) 

A=l 

where /j.^ are vectors of one of the fundamental representations (3 or 3) of SL{3). 

The particular form of integral representation ( [4.3.23 ) depends on particular screening 
insertions to the correlator ( [4.3.21|) . We will concentrate on the solutions which have no 
denominators. One of the reasons of such choice is that these solutions possess the most 
simple integrable structure, though the other ones can still be analyzed in the same manner. 
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The simplest solutions which have no denominators correspond to specific correlators 

^{CMMM,p-,i}(^^ = (n|e^(W) nexpQ^jO) (4.3.27) 

i 

when we take to be "neighbour" (not simple!) roots: (cKjCKj) = 1. In the case of 5*^(3) 
this corresponds, say, to insertions of only Qa^ and and gives formula (|4.2.30| ). Thus, we 



establish that the partition function ( 4.2.3(J| ) satisfies Pi^'^'^J'-constraint algebra, the generators 



of which can be read off from the formula ( [4.3. 19| ): 



wi'^ = 3Y.km-^ + 3Y.kt, gfgr+ (4.3.28) 

There are many possible generalizations of the above scheme, some of them, including 
" supersymmetric matrix models", can be found in |^ . 



Thus, we derived matrix models possessing PF-symmetry. The question we are going to 
discuss to complete this subsection is what symmetry does possess the standard multi-matrix 
model ([4.2.231) . Let us consider the simplest case of two- matrix model ([1.2.23[) , one of the 
potentials, say, V2 being truncated to be just cif|. Now we repeat the procedure of the 
derivation of the Virasoro constraints ( [4.3. 1[) . Let us do the following infinitesimal change of 



integration variables 

6H2 = Hi g > 

(4.3.29) 

5Hi = 3cHi {V({Hi) — H2) + "quantum corrections". 
This variation of variables induces the variation of potential Vi: 

SVi = -3c{V'fHl + Hl^^ + "quantum corrections" . (4.3.30) 

The first term in this expression gives rise to the "classical" part of the symmetry generators 
(analogous to above) and the second one produces the derivative d/dtg+i in the Ward 
identities. 



While the if2-component of the variation ( [4.3.29| ) does not change the integration measure 



dHidH2, this is not true for i72-component. The corresponding Jacobian is responsible for 
the "quantum" contributions to ([4.3.29| ) and ([4.3.30[) . 



Now we just write down the result for the constraints imposed on the partition function 
([4:2:301 ): 



dt 



^ +3cWil\{t}]zPit) 



= 0. (4.3.31) 

V2=cH^ 



where generators 



will = Y ktdk + Yktk E 



k,l>0 fc>o a+b=k+q~2 ^'^adh 

q(^q — X) Q 



(4.3.32) 



fc>0 a+b=q-l ^^^^^k+b-l a+b+c=q-2^^o.dtbdtc 2 dt 



1 
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produces some new constraint algebra. 

This W^^^ algebra arises in many different contexts and can be trivially generalized to 
any spin. The generators of this algebra have the same "classical part" as the standard 
W^^^ algebra (|4.3.28|) , but "quantum parts" are different. The set of VF-constraints is closed 
(as well Ward identities (|4.3.31| )) and is an algebra quadratic in generators, similar to the 
standard VT-algebra. Let us stress that considering the polynomial potential V2 of a degree 
K, one obtains the Ward identity, analogous (^4.3.31|) , but for W^^^-algehia. In particular, 
the generic potential leads to linear W^°^^-algehia. The detailed description of the properties 
of ly-algebra and constraints in two-matrix model (|4.2.23| ) can be found in ^ . 

Thus, we realized that the simple multi-matrix models (^4.2.23 ) satisfy unusual W^-constraints, 
while more complicated models ( [4. 2. 301 ) satisfy the standard W^- constraints. In the contin- 
uum limit this difference disappears, and both these types of models have the same (double 
scaling) continuum limit. Unfortunately, we do not know the efficient way of taking the con- 
tinuum limit in the multi-matrix models ( [4.2.23|) (certainly, I mean by this that we have no 
complete understanding of the renormalization of all operators and constraint algebra, i.e. of 
simultaneously all correlation functions). But CMMM turns out to be more suitable for the 
continuum limit procedure. This is one of their advantages. We are going to discuss these 
questions in the next subsection. 



4.4 Continuum limit 

4.4.1 Continuum limit of Toda chain equations 

In this subsection we duscuss the bridge between above considered discrete matrix models and 
these are in the double scahng limit. Let us stress that, due to singularity of this procedure, 
it is very difficult to take this continuum limit directly in the partition function (in particular, 
there should be some (infinite) renormalizations), but easier to do in some equations, where 
the porcedure is less singular (to be honest, there is a very simple way of doing the double 
scaling limit immediately in the partition function, but it requires some additional information 
and will be discussed in the next section). Indeed, we already mentioned that the crucial 
properties of integrability and imposed constraint algebra are preserved in the continuum 
limit. It implies that the most natural choice of equations to take the continuum limit would 
be either constraint algebra, or integrable equations. Now we work out both these ways for 
the simplest example of the Hermitian one matrix model (|4.2.1[ ). 

Let us start with integrable equations of Toda chain hierarchy. Our further discussion 
is sufficiently general and does not refer to the concrete matrix model r-function, i.e. the 
continuum limit conserves all the space of solutions to integrable equations, and the structure 
of integrability. From the other hand, we demonstrate below that it also conserves constraint 
algebra. It is specific of the double scaling limit, which can be invariantly determined just by 
its properties of not disturbing all crucial matrix model properties 

Thus, now we build the map from Toda hierarchy to KdV hierarchy by limiting procedure, 
and then we do analogous map from discrete Virasoro algebra ([4.3.3|) to that in the double 
scaling limit. Certainly, we consider here only the first equation of the hierarchy ( [4.2. 12[ ). For 
the sake of simplicity, we consider here only the case of even potential in ( |4.2.1| ). It implies 
that all odd times t2k+i = and all ipn are supposed to be independent of them. Thus, 
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p„ = = and the first non-trivial equation (which is really the second equation of Toda 
chain hierarchy) has a form (t = ^2): 

^ = -Rn{Rn+l - Rn-l). (4.4.1) 

This hierarchy is the reduction of Toda chain hierarchy and called Volterra hierarchy. 

Indeed, we could guess that the continuum limit of Toda chain hierarchy should be KdV 



hierarchy by looking at the Lax operator (|4.2.8|) . Indeed, the main idea of the continuum 



limit in matrix models is to change the discrete variable n by the continuous time variable 
(which will be the first time in the continuum hierarchy). Therefore, choosing n = xe we can 
get in the continuum limit: 



cont 



2 - e'[dl + (d^ix))]. (4.4.2) 

One can recognize in this operator just Lax operator of KdV hierarchy, the eigenvalues of 
Toda Lax operator being distributed around zero in the continuum limit (see 0]). Still for 



the equations the continuum limit procedure is not so simple. Indeed, doing naively, we get 



from (^4.4.11) : 

dR{x) 



-R{x){R{x + e) - R{x - e)) = -R{x')R'{x'). (4.4.3) 



dt 

provided by the replace x = ex'. This formula is called dispersionless KdV equation (or 
Bateman, or Hopf, or Khohlov-Zabolotskaya equation). This equation is too simple and 



corresponds to just naive "large n limit" of matrix models To obtain KdV equation, one 
should take care of non-linear term. The procedure allowing one to preserve such terms is 
called "double scaling limit" and is as follows. 

Imagine, that in continuum limit Rn tends to a constant Ro, and the function r{x) arises 
only as scaling approximation to this constant: R{x) = Ro{l + e^r{x)). Then the leading term 
at the r.h.s. of (g^) is eRR'{x) = -2eV(x)(l + 0{e\e')), and instead of ( ^X3| ) we would 
get: 

- = -26i?or'(x)((l + 0(6^e^)). (4.4.4) 
Then, by a simple change of variables^ 

Ti = x- 2eRot, Ts = e^Rot (4.4.6) 



it can be transformed into 



e-^0{e\e'), (4.4.7) 
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This change of variables is implied by the relation: 

3^ + 2^ - ( ^ + 26 i?o^ 1 + 2^ ^0^1 ^ - (4.4.5) 



26 



and terms at the r.h.s. also should be taken into account. Then we get: 



dr{x) 
dt 



-2eRo (^1 + eV(x))(r'(x) + ^eV'"(x) + 0(e^)) 
and, after the change of variables ( [4.4.6|) , 



(4.4.8) 



= -^r'"(Ti) - 2e^-Vr'(Ti) + 0{e\ e'). (4.4.9) 

One can see that the choice s = 2 is distinguished (a critical point) and at this point we get 
nothing but KdV equation: 

dr 1 d^r dr 

2rT^, (4.4.10) 



dTs 3 (9Tf (9Ti ' 

Thus, we obtained KdV equation. Certainly, one can get all KdV hierarchy in the same 
way. Moreover, one can investigate the continuum limit of the Toda (not Volterra) hierarchy, 
producing two equivalent KdV hierarchies (which correspond to flows in odd and even Toda 
times) - see |^0|, It reflects in some problems in the continuum limit of Virasoro algebra 
(see below). 

Thus, the lesson, which we can extract from this continuum limit procedure is that, to 
get the proper limit, one should do linear transformation of times and renormalization of the 
partition function (or i?„). Now we use these lessons in order to study the continuum limit 
of Virasoro algebra. 



4.4.2 * The continuum limit of the Virasoro algebra 

The procedure of taking the continuum limit of Virasoro algebra we demonstrate now is a bit 
tedious and too technical, as, in contrast to our previous consideration allows one to work 
with whole hierarchy simultaneously. Let me just describe a set of rules which leads to the 



correct answer, all details and hints can be found in 62 



It has been suggested in P^j that the square root of the partition function of the continuum 



limit of one-matrix model is subjected to the Virasoro constraints 

£conty^ = 0, q> -1, (4.4.11) 



where 



- + o) ^2fe+i^;^^ hC ^ 

(4.4.12) 

16 (16G) 
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are modes of the stress tensor 



r(.) = i:a*=W:-J3=i:-^. (4.4.13) 

Indeed, we demonstrate that these equations which reflect the Pi^^^J'-invariance of the 
partition function of the continuum model can be deduced from analogous constraints in 
Hermitian one-matrix model by taking the double-scaling continuum limit. The procedure is 
as follows. 

In order to obtain the relation between Ty-invariance of the discrete and continuum models 
one has to consider a reduction of model ( [4.2. 1| ) to the pure even potential t2fc+i = 0. Let us 



denote by the r^'^'^ the partition function of the reduced matrix model 

^"{^24 = / [dH] exp Tr ^ hkH^'' (4.4.14) 

k=0 

and consider the following change of the time variables 

9m — ^ ,N (4.4.15) 

k>m {k - my.T {m + 

where Qm = ^t2m and this expression can be used also for the zero discrete time qq = n that 
plays the role of the dimension of matrices in the one-matrix model. Derivatives with respect 
to t2k transform as 



dt 



2m+l 



where the auxiliary continuum times T2m+i are connected with "true" Kazakov continuum 
times T2m+i via 



^2fc+i = T2k+i + Q ,^ , 1 /o -^2(fc-i)+i, (4.4.17) 



k_ 

'k + l/2' 

and coincide with T2m+i in the double-scaling limit when a — 0. 

Let us rescale the partition function of the reduced one-matrix model by exponent of 
quadratic form of the auxiliary times T2m+i 



f = exp I -- A^kf2m+lf2k+l 1 (4.4.18) 



■m,k>0 



with 



r f A; + I) r fm + I) (_\k+m -k-m-l 

Au = ^ — - f4 4 19) 

212 (i) k\m\{k + m + l){k + m + 2)' v-- ; 

Then a direct though tedious calculation [B2l demonstrates that the relation 



28 



f 



g+1 

-q \ ^ 

2^ '-'q+l 



T red_red 



T 



red 



(4.4.20) 



is valid, where 



and 



L 



red 
2q 



k=0 



d 



2k 



dt 



+ E 



2(fc+g) 0<fc<q dt2kdt2(g-k) 



(4.4.21) 



fc>l 



k>0 



1 



2(A:-1)+1 

d 



16G'' 



fc>0 ^ 0'J2(fe+<7)+l 



9 



d 



Here = 



0<^-l 9T2k+l dT2(q-k-l)+l 

are binomial coefficients. 



n > 1. 



(4.4.22) 



These Virasoro generators differ from the Virasoro generators ( [4.4. 121 ) |T^ by terms 
which are singular in the limit a — > 0. At the same time U2p'^^'^'^ ^^e r.h.s. of ( [4.4. 20|) do 



not need to vanish, since 



= L2pT 



T red red 
^2p ^ 



E 



dt2i+idt2 



(4.4.23) 



i2k + l—0 



It was shown in |^| that these two origins of difference between ( 4.4. 12| ) and ( [4.4.22| ) actually 
cancel each other, provided eq. ( [4.4.20|) is rewritten in terms of the square root \^ rather than 
f itself: 



/'cont q+1 T 

p=0 



t2fc+l=0 



;i + 0(a)) 



(4.4.24) 



The proof of this cancelation, as given in |£2[, is not too much simple and makes use of 
integrable equations for r. 



4.4.3 * Invariant formulation of the limiting procedure 

After we demonstrated the continuum limit procedure for the Hermitian one matrix case, 
the next question to be addressed to is how to generalize this limiting procedure to multi- 
matrix case. Indeed, in the case of conventional multi-matrix models this procedure is still 
unknown. But for CMMM it can be trivially done, what is one of the advantages of CMMM we 
mentioned above. To do this, let us demonstrate a more economical way to define the change 
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of the time- variables t — > T (also discussed in [§2]), implied by the scalar field formalism. 
The Kazakov change of the time variables ( [4.4.15| , ^.4.16|) can be deduced from the following 



prescription. Let us consider the free scalar field with periodic boundary conditions 

d^{u) = E g.u^'^'' + E j^u-''~\ (4.4.25) 

fc>0 k>l '^''2fc 

and analogous scalar field with antiperiodic boundary conditions: 

dm = Y.i(k + l) T2k^i^'-^ + T^^-'-O ■ (4.4.26) 
fc>oV^ dT2k+i J 



Then the equation 



ld<i>{z)f = a-^d^iuy""^, = l + az (4.4.27) 



generates the correct transformation rules (|4.4.15| ), ([4.4.16|) and gives rise to the expression 



for Afcm ( [4.4.191 ). Taking the square of the both sides of the identity ( [4.4.27| ), 



^T{z)f 

T 



T 



red 



T(u)t 



red 



(4.4.28) 



one can obtain after simple calculations that the same relation ( |4.4.2(]| ) is valid. 

Now it is evident how to generalize this procedure to general multi-matrix models using 
the formalism of scalar fields with Zp-twisted boundary conditions. It can be found in details 
in the references fff^, |49 



4.5 Fermionic representations and forced hierarchies 
4.5.1 Determinant representations of r-functions 

Beginning with this subsection we are going to spend more time discussing integrability in 
matrix models, which was practically out of our consideration before. In particular, now we 
are going to discuss the notion of r-function, main definitions and the determinant representa- 
tions, which we were talking about in the second subsection. We mainly follow the references 

1631, m ra. 



We restrict ourselves to consideration of ordinary KP and Toda-lattice r-functions, which 
can be constructed as correlators in the theory of free fermions ipjijj {b, c-system of spin 1/2). 
The basic quantity is the ratio of fermionic correlators. 



{N 


e^Ge 


H 


N) 


{N 


G 


N) 



rNit,t\G) = ' ^, ' ^ (4.5.1) 



in the theory of free 2-dimensional fermionic fields ip{z), ip*{z) with the action / ip*d'ip (see 
the Appendix Al for the definitions and notations). In the definition ( [4.5. 1|) 
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H=J2tkJk; H = J2hJ^k, (4.5.2) 

fc>0 A;>0 

and the currents are defined to be 

j{z) = r{zm^); (4.5.3) 

where 

^(z) = 5:^„^" dz'/\ 

^ , (4.5.4) 

z 

This expression for r-function is rather famous. Let us point out that this r-function is the r- 
function of the most general one-component hierarchy which is Toda lattice hierarchy. Putting 
all negative times to be zero as well as discrete index n, one obtains usual KP hierarchy. An 
element 

G=:exp{Y,Gmn^M: (4.5.5) 

m,n 

is an element of the group GL(oo) realized in the infinite dimensional Grassmannian. The 
normal ordering should be understood here with respect to any \k) vacuum, where vacuum 
states are defined by conditions 

^^\k) = m<k , ^*Jk) = m>k. (4.5.6) 

Any particular solution of the hierarchy depends only on the choice of the element G (or, 
equivalently it can be uniquely described by the matrix Qkm)- Thus, we can reformulate the 
sense of constraint algebra in matrix models: it just fixes the element of the Grassmannian, 
i.e. one can reformulate constraint algebra in terms of equations on the element of the 



Grassmannian. This approach was advocated in |6J]. 

From commutation relations for the fermionic modes, one can conclude that any element 
in the form ( |4.5.5| ) rotates the fermionic modes as follows 



G^kG-' = ^,R,u , GruG-' = r.Rkh (4.5.7) 



where the matrix Rj^ can be expressed through Qj^ (see [^). We will see below that the 
general r-function (|4.5.1| ) can be expressed in the determinant form with explicit dependence 



of Rjk ■ In order to calculate r-function we need some more notations. Using commutation 

relations for the fermionic modes, one can obtain the evolution of ^"^^ ^*(^) i^i times 
{tk}, {ik} in the form 

i){z, t) = e^(*V(^)e"''^*^ = e«(*'^V(^) , (4.5.8) 

t) = e^(*V*(^)e-''^*) = e-«(*'^V*(^) ; (4.5.9) 
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Indeed, it is rather crucial point and we will return to possible choices of the normal orderings. 
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-Hit) 



ij{z) 



where 



i(t,z) = Y.tuz^ . 



k=l 



Now let us define Shur polynomials by the formula: 



(4.5.10) 
(4.5.11) 

(4.5.12) 



(4.5.13) 



k>0 k>0 

Using this definition, from eqs. ( |4.5.8| )-( [4.5.11| ) one can easily obtain the evolution of the 
fermionic modes: 



m=0 

oo 

rkit) ^ e^w^^e^w = E rk^M-t); 

m=0 

_ oo 

Mi) = e^^V-fce-^® = J2 ^k+mPrnii); 

m=0 
_ oo 

rkit) ^ e^(*")^:e-^® = Y: rk-mPU-i)- 

m=0 



(4.5.14) 
(4.5.15) 
(4.5.16) 
(4.5.17) 



It is useful to introduce the totally occupied state | — oo) which satisfies the requirements 



- oo) = , i G z . 
Then any shifted vacuum can be generated from this state as follows: 



(4.5.18) 



\n) = ip„-iip„. 



oo) 



(4.5.19) 



Note that the action of an any element G of the Clifford group (and, as consequence, the 
action of e~^®) on | — oo) is very simple: G\ — oo) ~ | — oo), so using ([4.5. 15|) and ( [4.5. 16|) 
one can obtain from eq.( |4.5.lD : 



r„(t,t) = (-oo|...<_2(-t)V^:_i(-t)G^„_i(-t)^„-2(-t)...| - oo) 
~ det[(-oo|^;(-t)G'^,(-t)G'-i| - oo)] . 
Using ( |4.5.7| ) it is easy to see that 



(4.5.20) 



Gr{-i)G-' =J2Pm{-i)rRk,j- 



(4.5.21) 



m,k 



and the "explicit" solution of the two-dimensional Toda lattice has the determinant represen- 
tation: 
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Tn{t,i) ~ det Hi 

+n,j+n 



i,j<0 



(4.5.22) 



where 



Hij{t,i) — RkmPk~iit)P„ 



(4.5.23) 



The ordinary solutions to KP hierarchy correspond to the case when the whole evolution 
depends only of positive times {tfc}; negative times {tk} serve as parameters which parameter- 
ize the family of points in Grassmannian and can be absorbed by re-definition of the matrix 
Rkm- Then r-function of (modified) KP hierarchy has the form 



Tn{t) = (n|e^(*)G'(t)|n) ~ det 
where G{t) = Ge"^® and 



Rk,j+n{t) Pk-i-n{t) 

. k 



(4.5.24) 



i,j<0 



Rkjiy) = ^RkmPm~j{-y) ■ 



(4.5.25) 



4.5.2 Determinant representations for matrix model hierarchies 

Now let us briefly discuss the formulas of the second subsection within the developed frame- 
work. Indeed, we already have the general determinant representation (|4.5.22| ) and would have 
to reduce it to concrete subhierarchies. The life is, however, not so simple as all the determi- 
nants we know from matrix models are finite determinants, in contrast to ( [4.5.22|) . This is an 
essential problem. Indeed, this is due to the specific boundary condition like ( [4.2. 15| ) imposed 
on matrix model r-functions. The hierarchies with this boundary condition are called forced 
hierarchy and studied in details in [^2|, The main result of this study is that these 
hierarchies are corresponded by the singular elements of the Grassmannian, which are induced 
by the fermion operators ( [4.5.4|) , but with only positive (or only negative) modes taken into 
account. It implies the consideration in ( 4.5.5|) of semz-infinite (instead of infinite) matrix 
of the group GL(oo). In its turn, it leads to the finite (instead of semi-infinite) determinant 



Indeed, there is the problem of continuing the boundary condition ( 4.2.15| ) to negative 
n. Say, the Toda chain integrable equations can not fix it (it can be trivially observed from 
Hirota equation - see ( |4.2.14| )). However, switching on the negative (Toda) times allows one 
to fix this continuation. And if one wishes to have the Toda chain as a proper reduction from 
Toda lattice hierarchy, it is necessary to put 







(4.5.26) 



(the case of r_„ = r„, for example, corresponds to CKP system - see |^). To do this, one 



should multiply the element of the Grassmannian by the projector onto positive n: 



exp[^ ^jjiijj* 

i<0 



(4.5.27) 
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with the properties: 

P^\n) = 9{n)\n), 
P+r-k = ^-kP+ = , A; > ; (4.5.28) 

The point, however, is that after this is done, the defining properties of the original (non- 
forced) system are the same as those of the forced one under the condition (|4.5.26| ) (i.e. 
after throughing away the half of fermionic modes and inserting the projector). Therefore, 
we discuss now the Toda lattice hierarchy and its Toda chain reduction with no restricting 
ourselves to the forced case. 



Let us note that the entries ( |4.5.23| ) of the determinant ( [4.5.22 ) are not arbitrary ones, 
but are conditioned to satisfy the properties: 



dHij/dtp = Hij^p, j > p>0, 
dHij/dtp = -Hi_pj, i> p> 0, 

which are the automatic consequences of the corresponding property of Shur polynomials: 



(4.5.29) 
(4.5.30) 



dPk/ dtp — Pk-p 



(4.5.31) 



following immediately from their definition. From the properties ( [4. 5. 291 ) and ( [4.5.30| ) one 
obtains directly our conditions (|4.2.19|) and ([4.2.21|) and their analogs for the multi-matrix 
case ( |4.2.27| ) (in the latter case, one should change the sign in the last potential in ( [4.2.23| ), 
which is fixed by the sign in the Hamiltonians in ik). Now let us discuss the condition ( [4.2.20|) , 
i.e. the reduction of general Toda lattice hierarchy to Toda chain. It can be easily formulated 
as condition to the element G of the Grassmannian: 



[Jk + Jk,G] = 0, 



which is equivalent to constraint 



[A + A-\R] = 0., 
where A is shift matrix Ajj = In this case, 

(^Q~^kJk — ^ — tkJk^ — tkJk(^ ^^kJk 

and r-function depends (up to the trivial factor) only on times {t^ —1^}: 
The reduction ( [4.5.33|) has an important solution^ 



(4.5.32) 



(4.5.33) 



(4.5.34) 



(4.5.35) 



^^Generally the solutions Rnk = Rn+k and Rnk = Rn-k are different, but for forced hierarchy, when r- 
function is the determinant of finite matrix, these two solutions are equivalent due to possibility to reflect 
matrix with respect to vertical axis without changing the determinant. 
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Rkm — Rk+m ■ (4.5.36) 



In this case the matrix Hij defined by eq. ([4.5.23|) evidently satisfies the relations Hij = Hi^j 
and 

((9(^ + dt^)Hi^j = for any k < n — i , k < n — j (4.5.37) 
due to the properties ( [4.5.29|) and (|4.5.30| ). The property ( [4.5.37|) certainly does not imply 



that the corresponding r-function depends only on difference of times because of restriction of 
values of k, but it restores correct dependence of times with taking into account of exponential 

in {w;^. 



Thus, we have proved that, for the Toda chain hierarchy, the matrix Hij satisfies the 
condition 

[H , A + A-'^] = 0, (4.5.38) 

which leads to r-function of Toda chain hierarchy (properly rescaled by exponential of bilinear 
form of times) which depends only on the difference of positive and negative times tp — tp, but 
not on their sum (one can consider this as defining property of Toda chain hierarchy). Let 
us remark that only one possible solution to constraint ( |4.5.38| ) is matrix Hij = Hi^j, but 



it should not be generally independent of the difference of times tp — tp. Certainly, one can 
through out negative times as the final answer for complete objects like r-function should be 
really independent of this difference. 

Anyway, we obtained for the r-function of Toda chain hierarchy the determinant repre- 
sentation with necessary conditions ( [4.2.19| )-( |4.2.21| ), i.e. the matrix Hij just corresponds to 
moment matrix Cij and we finally reproduce the result ( [4.2.18|) . 



4.5.3 * Integrable structure of CMMM - multi-component hierarchy 

Now we are going to say some words on integrable structure which arises in CMMM. In fact, 
we already wrote down the determinant representation ( [4.2.31|) for CMMM. Now there are 



two things which can be established. The first one is to obtain the integrable equations. It 
can be done immediately from the equation (^4.2.31|) and leads to clear generalization of the 
Hirota equation ( ^.2.14| ): 



tn M dtidti ' V / V 



The second point is just the properties of the r-function corresponding to CMMM within 
the framework described above. The statement is that the partition function of CMMM is a 
r-function of a multi-component ((p+l)-component) Kadomtsev-Petviashvili hierarchy which 
obeys the constraint 

X:9M'=V(^+^)({t}) = 0, n = l,2,..., (4.5.40) 

k=l 
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where p is the number of matrices in CMMM. Let us demonstrate how it works in our usual 
simplest example of p = 1. It should be 2-component hierarchy. 

The r-function of 2-component KP hierarchy is by definition the correlator 

4'm(^, i) = {N, M|e^(*'*"^G'|iV + M, 0) (4.5.41) 

where 

Hit,t) = Y.it,4'^+h4'^) (4.5.42) 

fc>0 

J»(^) = J2 Jk^z-'^'' = ■.iJ^'\z)ij^^^*{zy. (4.5.43) 

^« (2)^0> (z') = + ... . (4.5.44) 

All these formulas are trivially generalized to {p + l)-component case by introducing p + 1 
fermions. 



Now we are going to demonstrate that ( [4.3. 16| ) is equivalent to ( |4.5.41| ) for certain G for 



which ( [4.5.411) depends only on the differences tk — t^. To do this we have to make use of the 



free-fermion representation of S'L(2)fc=i Kac- Moody algebra: 

Jo = ^(^(l)^(^)* - ^(2)^(2)*) ^ 1 (J(l) _ J(2)) 

J+ = ^(2)^(1)* j_ ^ ^(1)^(2)* (4.5.45) 

Now let us take G to be the following exponent of a quadratic form (indeed, we know that 
this is the correct fermionic form for exponential of screening operators) 

G = : exp (^J V^^V^^^*) : (4.5.46) 

Now we bosonize the fermions 

jW = d(f)u J^^^ = d(f)2 (4.5.47) 

and compute the correlator 

- M<^' (g^'^^^'^ + ''■^^'^^) (/ ^V^^'V^*:)'" |0) = 

= ^(iV, -N\ exp (^f[V^{z)j('\z) + V2{z)J^'\z)]^ (| : exp(0i - 0^) |0), 

where V(i,2) = J^ki^k, ik)z^. Introducing the linear combinations -\/20 = 0i— 02, V^<P = 0i+02 
we finally get 
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x(iV|exp (l=j^[V,{z)-V2{z)]d^{z) \ (I :expV24>-y \0)=rg\t-i) 
since the first correlator is in fact independent of t and t. Thus, we proved that the r-function 



( 4.5.41 ) indeed depends only on the difference of two sets of times {tk — tk}- So, we obtained 



here a particular case of the 2-component KP hierarchy (|4.5.41| ) and 

(i) requiring the elements of Grassmannian to be of the form ( [4.5.46|) we actually performed 



a reduction to the 1-component cas^; 

(ii) we proved in ( |4.5.48| ) that this is an AKNS-type reduction for the r-function ( [4.5.47|) 



Mm- 

The above simple example already contains all the basic features of at least all the Ap cases. 
Indeed, the reduction ( [4.5.48|) is nothing but S'L(2)-reduction of a generic GL{2) situation. In 
other words, the diagonal U{1) GL(2)-current J = |( J*-^^ + J^^-') = -^dcp decouples. This is an 
invariant statement which can be easily generalized to higher p cases. Indeed, the screening 
operators for the p-matrix case are the integrals of SL{p + l)i Kac-Moody currents (not 
GL{p + 1) ones) and, thus, the r-function ([4.3.21|) - ([4.3.22|) does not depend on x^^}, 
i.e. we obtain the constraint ( |4.5.40| ). 



As to the Toda-like representation of CMM, in the simplest S'L(2)-case the result should 
be equivalent to the Toda chain hierarchy. In the fermionic language this connection is 
established by the following substitution in the element of the Grassmannian 



'^^^\z) ^ipiz), iIj^^\z) ^ ipi-) (4.5.49) 

and the same for '?/'*'s. This is a reflection of the fact that Toda system is described by the two 
marked points (say, and oo) and corresponds to two glued discs, so it can be also described by 
two different fermions. This might lead to a general phenomenon, when any multi-component 
solution to CMM is actually related to (some reduction) of a multi-component Toda lattice. 



5 GKM approach to matrix models 



5.1 What are the matrix models in the double scaling limit 



In this section we are going to discuss the models which arise after the continuum (double 
scaling) limit is taken. Indeed, now we are merely formulate a set of conditions, which uniquely 
defines the corresponding partition function. 



Note that the idea to preserve both indices in ( 4.5.41 ) leads immediately to additional insertions either 
of ^p'^^^* or ip'^^^ to the right vacuum |0), so that it is no longer annihilated at least by the r_i Virasoro 



generator, or in other words this ruins the string equation. Thus only the particular reduction (4.5.46) seems 
to be consistent with string equation. This choice of indices just corresponds to that considered originally in 
f|6l. 
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These conditions are of two kinds. The first one describes the integrable properties. That 
is, the square root of the partition function of the Hermitean i^'-matrix model in the continuum 
hmit is a r-function of the iT-reduced KP hierarchy. Say, one-matrix model corresponds to 
KdV hierarchy, two-matrix model - to Boussinesq hierarchy and etc. 

The second condition imposed on the partition function is so-called string equation. The 
role of this equation is to fix the r-function which corresponds to the matrix model. Put it 
differently, this equation gives the point of the Grassmannian which corresponds to the matrix 
model. The equation is 

^^^-Vr = 0, (5.1.1) 

where 



a,b>0 

Indeed, hereafter we call string equation a bit more general equation 



Cifr = 0, 



(5.1.3) 



which is also correct. 



It can be demonstrated that above mentioned conditions imply that the r-function 
satisfies a whole set of constraints. More concretely, i^-reduced r-function satisfies the set of 
W-algebra constraints: 



W«r^^> = 0, 



(5.1.4) 



l^-operators here are the standard generators of Fateev-Lukyanov-Zamolodchikov 14^-algebra 
5B| expressed in the manifest terms of times, i.e. with creation and annihilation operators 
being T„ and d/dTn respectively and acting on the spce of function of times |]32[| . 

There is also inverse conjecture which asserts that the constraints (|5.1.4| ) have a unique 
solution, i.e. they select out the r-function of i^-reduced KP hierarchy, which satisfies the 
string equation (|5.1.3| ). This conjecture is still not proved, though there are strong arguments 
in favor of its correctness 



Let us note now that it is not so simple to prove the properties described above, and, 
in part, they are merely an invariant definition of the double sclaing continuum limit. Some 
general argument have been used to derive these properties |T^, ^], but there are still only 
two rigid calculations. 

The first one describes the transition from discrete to continuum W- and Virasoro con- 
straints and has been done for ii" = 2 case of the standard matrix models p2] and for the 
general case of the conformal multi-matrix models ||4^, Another calculation makes use 
of the manifest matrix integral representation of the partition function |^5|, which will be 
described in the next subsections. 
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5.2 Generalized Kontsevich Model (GKM) 



Thus far we have described the set of constraints which uniquely defines the partition function 
of matrix models in the double scaling limit. Now we are going to present the manifest solution 
of these constraints. Indeed, for the simplest case of one-matrix model it was derived from 
Witten topological theory (see section 3) by Kontsevich |6^. We propose proper generalization 
of his result, demonstrating that our solution satisfies all necessary constraints (in contrast 
to the approach where there has been manifestly calculated the integral ( p.l[ )). 



For the lack of space, we will briefly describe only the structures and appealing properties 
of the partition function, which we call Generalized Kontsevich's Model (GKM) |]7D|. Unfor- 
tunately, we can not explain here the deep connections of GKM with the theory of integrable 
systems and develop possibile generalizations, valuable in different applications. In particular, 
the connection with topological N = 2 Landau- Ginzburg theories is out of the scope of the 
present review. We refer to the conclusion for a short review of the existing literature. 

The partition function of the GKM is defined by the following integral over N x N Her- 
mitean matrix: 

Z^M] - J^trn^ . (5.2.1) 

where 

f/(M, Y) = Tr[V{M + Y)- V(M) - V '{M)Y] (5.2.2) 

and 

U2(M, Y) = lim ^U(M, eY) , (5.2.3) 

is an y^-term in U. M is also a Hermitean N x N matrix with eigenvalues {/Xj}, V(/i) is 
arbitrary analytic function. 

All crucial properties of GKM follows from its definition. That is, its partition function 
may be considered as a functional of two different variables: potential V(/i) and the infinite- 
dimensional Hermitean matrix M with eigenvalues {/ij}. Partition function is an A^- 
independent KP r-function, considered as a function of time-variables T„ = ^TrM~^ and the 
point of Grassmannian is specified by the choice of potential. The A^-dependence enters only 
through the argument M : we return to finite-dimensional matrices if only A^ eigenvalues of 
M are finite. In this sense the "continuum" limit of A^ — > cxd is smooth. 

The GKM is associated with a subset of Grassmannian, specified by additional 
constraint. For particularly adjusted potentials V(/i) = const- fi^~^^, the corresponding points 
in Grassmannian lies in the subvarieties, associated with i^'-reductions of KP-hierarchy, Z^^^ 
becomes independent of all the time-variables T^n, and the £_i-constraint implies the whole 
tower of H^i^-algebra constraints on the reduced r-function. These properties are exactly the 
same as suggested for double scaling limit of the K — 1-matrix model, and in fact there is an 
identification 



ZL^' = \/rif-^>. (5.2.4) 
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All this means that GKM provides an interpolation between double-scaling continuum 
limits of all multimatrix models and thus between all string models with c < 1. Moreover, 
this is a reasonable interpolation, because both integrable and "string-equation" structures 
are preserved. This is why we advertise GKM as a plausible (on-shell) prototype of a unified 
theory of 2d gravity. 

In the next subsections we briefly comment all properties described here. 

5.3 Main integrable properties of GKM 
5.3.1 Integrable structure 

After the shift of variables X = Y + M and integration over angular components of X, 
zj^\M] acquires the form of 

zPm = (5.3.1) 

where A(M) = — Hj) is the Van-der-Monde determinant, and functions 

|.,(^) = [V "(/x)]i/2e^('^)-'^^ '^^^ f e-^(")+"^'(^)x*tix (5.3.2) 



The only assumption necessary for the derivation of (4) from (1) is the possibility to represent 
the potential V(/i) as a formal series in positive integer powers of /i. 

Formula ( |5.3.1|) with arbitrary entries 0j(/i) having the asymptotics 

(Pifir^ fx'-\l + 0{fx)) (5.3.3) 
is characteristic for generic KP r-function T^(Tn) in Miwa's coordinates 

T„ = -TrM-"" , n > 1 (5.3.4) 

n 

and the point G of Grassmannian is defined by potential V through the set of basis vectors 
{0i(/i)}. The most immediately this set is connected with the Baker- Akhiezer (BA) function 

v^±(.|T,)=eS^-^^i^^^^ (5.3.5) 

through the relation: 



^+(/i|rfe 



Thus, we obtain 



ZW[Af] =r^^}(T,). (5.3.7) 
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The case of finite N in this formalism is distinguished by the condition that only of the 
parameters {ni} are finite. In order to take the limit — > oo in the GKM ( ^.2.1|) it is enough 
to bring all the fi'^s from infinity. In this sense this a smooth limit, in contrast to the singular 
conventional double-scaling limit, which one needs to take in ordinary (multi)matrix models. 

Let us check that this limit is really smooth. To do this, one should change the number 



of Miwa variables + 1 — bringing fiN+i to oo in the determinant ( |5.3.1| ). Using the 



asymptotics ( 5.3.3| ), it is trivially to prove that it results in A^ x A^ determinant of the same 
type. This means that A^ enters only as the number of Miwa variables, but not as an explicit 
parameter. 



5.3.2 Reductions 

The integral jF'f^^[A], A = V '(M), in the numerator of ( |5.2.1| ) satisfies the Ward identity 



Tr I e(A) 



V'( 



d 



OA 



-A 



tr 



{V} 



N 







(5.3.8) 



(as result of invariance under any shift of integration variable X — >■ X + e(M)). Let us denote 
the integral in ( |5.3.2| ) through Fj(V'(/i)). If potential V(/u) is restricted to be a polynomial of 
degree i^T + 1, this identity implies that the functions -Fj(A) obey additional relations: 



(A) = A" ■ F„(A) 



m+Kn—l 
i=l 



(5.3.9) 



Since the sum at the r.h.s. does not contribute to determinant ( |5.3.1| ), we can say that all the 
functions Fn are expressed through the first K functions Fi...Fk by multiplication by powers 
of A = V '(/i). Such situation (when the basis vectors 0j, defining the point of Grassmannian 
are proportional to the first K ones) corresponds to reduction of KP-hierarchy. This reduction 
depends on the form of V '{jj) and in the case of V^ji) = VxifJ') = const ■ /x^"*"^ coincides with 
the well-known i^-reduction of the KP-hierarchy (KdV as K = 2, Boussinesq as = 3 etc.). 
Thus in such cases partition function of GKM becomes r'>^^J'-function of the corresponding 
hierarchy. Generic t^^^ possesses an important property: it is almost independent of all 
time- variables T„x- To be exact, 

d logT^^^/dTnK = a„ = const. (5.3.10) 

In variance with generic t^^^\ the partition function Z^^^ of GKM obeys this condition with 
all ttn = (see proof in |7D[|). 



5.4 String equation 
5.4.1 £ i-constraint 

The set of function {$j(/i)} in (4) is, however, not arbitrary. They are all expressed through 
a single function — potential V(/x), —and are in fact recurrently related: 
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F,(A) = {d/dXy-'Fi{X). 
This relation is enough to prove that 



d 



{V} 



d 

-Tr M + Tr—— log det Fi{Xj] 

uh-tr 



whenever potential V(/i) grows faster than /x as /x — oo. 
Thus, Z^^^ satisfies a simple identity: 



1 



^{V} y{V} 



where operator C^i is defined to be 



— log + TrM - Tr-— log det F,(A,) = 

OI 1 OAtr 



(5.4.11) 



(5.4.12) 



(5.4.13) 



1 



d 



2 tr ^ "(/^^)'^' "(/^.o 



Hi - fij 



dTi 



(5.4.14) 



(the items with i = j are included into the sum). From eqs.( |5.4.TB| ),( p.4.T^ it follows that 
partition function of GKM satisfies the constraint 



(5.4.15) 



The reason why the operator ( ^.4.14 ) is denoted by is clear from the following argu- 
ment. If V = Vk, the generic expression ( |5.4.14| ) for the £_i-operator turns into 



^ n>K a+b=K 



(5.4.16) 



a,b>0 



The last item at the r.h.s. may be eliminated by the shift of time-variables: 

K 

n 

Only expressed in terms of these T's the constraint ( [5.4.14 ) acquires the form of 

^ n>K a + b=K 

y n^OmodK a,6>0 / 



(5.4.17) 



(5.4.18) 



coinciding with expression ( |5.1.3|) . 
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5.4.2 Universal string equation 

Generalization of the string equation to the case of arbitrary potential 
may be transformed to the following form 



. 



d'^ log r 



n>- 



1 dT.dTn 



u , 



where 



(5.4.19) 



(5.4.20) 



u 



% = Tr 



1 



1 



V"{M) M" 



+1 



d"^ log r d log T _ d log t _ „ 



on 



dT_i 



(5.4.21) 



Using BA function (|5.3.5| ), one can rewrite string equation (|5.4.15|) in the form of bilinear 
relation 



(5.4.22) 



5.4.3 W-constraints 

According to subsection 5.1 the constraint 

implies the entire tower of W-constraints 



(5.4.23) 



0, k = 2,3, ...,K: n> 1-k 



(5.4.24) 

imposed on t^^\ Here is the n—th harmonics of the p—th generator of Zamolodchikov's 
H^x-algebra (the proper notation would be W^p^^^\ but it is a bit too complicated). There 
is a Virasoro Lie sub-algebra, generated by Wj^^ = Cl^^\ and the particular is just the 
operator (|5.4.16|) . This is the true origin of our notation in the generic situation (where 



the entire Virasoro subalgebra of W^o was not explicitly specified). 

Besides being a corollary of ( [5. 4. 231 ), the constraints ( |5.4.24| ) can be directly deduced from 
the Ward identity ( |5.3.8| ). For the case of = 2 (which is original Kontsevich's model []69| ) 
this derivation was given in (see also|]72 



j, for alternative proofs). Unfortunately, 
for K > 3 the direct corollary of ( ^.3.8|) is not just ( 5.4.24|) , but peculiar linear combinations 
of these constraints, e.g. for K = 3 they look like 
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WlJzW = 0, n > -2; 



fc>i 



9 



= 0, 



Cirp ' '36—3 ( oo 

a+6=3n "^"^ 3^+2 



a,b>0,n> -2; 
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E (3A: - 2)T3,_2W^fV3n + E 

_fe>l+6„+3,o a+fe=3n^^3a+l 



)A;{2) I v{3} _ Q 



(5.4.25) 



a,b > 0, n > —3. 

For identification of ( |5.4.25| ) witli ( |5.4.24| ) one can argue tliat botfi sets of constraints possess 
unique, and tlius coinciding, solutions. Really, only the case of i^T = 3 has been investigated 
in this case in the paper 



6 Conclusion 

In the present review we have discussed in some details only Liouville (physical) gravity and 
discrete matrix models. We almost ignored topological theories, and described GKM approach 
very briefly, though they deserve separate reviews. 

Indeed, this is the GKM approach which allows one to connect different branches and 
viewpoints. On one hand, for the particular case of the cubic potential {K = 2) GKM integral 
can be immediately derived from the Witten topological theory |jl8|]. On the other hand, it 



is possible to establish connection between GKM and N = 2 Landau- Ginzburg topological 



gravity [^6|, As topological theories of both these types are appropriate models of 2d 
gravity, the connections of these with GKM are of great importance. 

Let us note that there are many generalizations of the original GKM integral ( |5.2.1| ). Say, 
one can include "negative" and "zero" time variables to obtainToda lattice r-functions |^^. It 
allows to incorporate in GKM treatment Hermitean one-matrix model and to take the double 
scaling limit immediately in the partition function . 

Another possible generalization is to consider polynomials of as the GKM potential. 
This case has much to do with unitary matrix integrals . 



An absolutely different issue is the existence of different phases in GKM. The new one can 
be defined by an other Miwa transformation of times: 

fk = ^TrM^ (6.0.1) 
k 

In fact, it rather determines a different aymptotics of GKM integral than a new phase. It 
acquires the sense of phase only when GKM integral is applied to the description of QCD 



A separate issue is the connection between GKM integral and character expansion of r- 
functions [^]. This point is rather important in applications both to unitary matrix models 
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and to 2d QCD 17|, |68[.g 

Another limitation of this review is that we restricted ourselves to Hermitian only matrix 
models and only to one-cut solution in the continuum limit. In fact, the properties, say, of 
unitary matrix model fSO, l78|, |81], ET[ is understood well too. Its integrable properties are 



described at discrete level |^2[ as well as in the continuum limit | 
algebra was also investigated, again both at the discrete [B^, pT| 



^, |8J, |79[ . Proper Virasoro 
and continuum levels. 



Different reductions were discussed |£5[ and GKM framework proposed |7^, ^ . 

Nevertheless, these models have no clear physical implications for 2c? gravity, but rather in 
2d QCD. This model is also not included into physically essential general framework, which 
would be multi-component hierarchies. These have (presumably) to include at equal footing 
DKP hierarchy and some other reductions having to do with different 2d gravity theories, 
including those with c = 1. These models should be associated with non-trivial critical points 
like multi-cut solutions iSBI . 



Thus, unitary model, on one hand, requires a separate review and, on the other hand, is 
out of our main line. This is why we avoided to discuss them here. 

This discussion might explain that it is impossible to put all the material on the small 
room of this review. The above-mentioned gaps will be filled in the second part of the review, 
devoted exclusively to matrix models in external fields, which have much to do both with 2d 
gravity (see section 5 above) and 2d gauge theories. 
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